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Abstract 

In this paper, solutions to the Ernst equation are investigated that de- 
pend on two real analytic functions defined on the interval [0, 1]. These 
solutions are introduced by a suitable limiting process of Backlund trans- 
formations applied to seed solutions of the Weyl class. 
It turns out that this class of solutions contains the general relativistic 
gravitational field of an arbitrary differentially rotating disk of dust, for 
which a continuous transition to some Newtonian disk exists. It will 
be shown how for given boundary conditions (i. e. proper surface mass 
density or angular velocity of the disk) the gravitational field can be 
approximated in terms of the above solutions. Furthermore, particu- 
lar examples will be discussed, including disks with a realistic profile 
for the angular velocity and more exotic disks possessing two spatially 
separated ergoregions. 
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1 Introduction 



Differentially rotating disks of dust have already been studied by Ansorg and Meinel 
|jj . They considered the class of hyperelliptic solutions to the Ernst equation intro- 
duced by Meinel and Neugebauer ||, see also j3|-||. These hyperelliptic solutions 
depend on a number of complex parameters and a real potential function. Ansorg 
and Meinel concentrated on the case in which one complex parameter can be pre- 
scribed. They determined the real potential function in order to satisfy a particular 
boundary condition valid for all disks of dust. To generate their solutions, they 
used Neugebauer's and Meinel's rigorous solution jjj || ||] to the boundary value 
problem of a rigidly rotating disk of dust which also belongs to the hyperelliptic 
class. 

A subclass of Ansorg's and Meinel's solutions is made up of Backlund transforms 
of seed solutions of the Weyl classQ. Solutions of this type are of particular interest 
since their mathematical structure is much simpler than that of the more general 
hyperelliptic solutions. 

With this in mind, the following questions arise: 

• Is it possible to find solutions corresponding to more general differentially 
rotating disks of dust by increasing the number of prescribed complex param- 
eters? 

• If so, is there a rapidly converging method for approximating arbitrary dif- 
ferentially rotating disks of dust with given boundary conditions (i. e. proper 
mass density or angular velocity)? 

• Is it perhaps possible to construct such a method by restriction to the much 
simpler solutions of the Backlund type? 

To answer these questions, the paper is organized as follows. In the first section 
the metric tensor, Ernst equation, and boundary conditions are introduced and the 
class of solutions of the Backlund type is represented. As will be discussed in the 
second section, the properties of these solutions can be used to obtain more general 
solutions by a suitable limiting process. Since these more general solutions depend 
on two real analytic functions defined on the interval [0,1], a rapidly converging 
numerical scheme to satisfy arbitrary boundary conditions for disks of dust can be 
created. This is depicted in the third section. Finally, the fourth section contains 
particular examples of differentially rotating disks of dust, including disks with 
a realistic profile for the angular velocity and more exotic disks possessing two 
spatially separated ergoregions. 

In what follows, units are used in which the velocity of light as well as Newton's 
constant of gravitation are equal to 1. 



1.1 Metric Tensor, Ernst equation, and boundary conditions 

The metric tensor for axisymmetric stationary and asymptotically flat space-times 
reads as follows in Weyl-Papapetrou-coordinates (p, £, ip, t): 

ds 2 = e~ 2U [e 2k (dp 2 +d( 2 )+ p 2 d(p 2 } - e 2U \dt + a dip) 2 . 

For this line element, the vacuum field equations are equivalent to a single complex 



equation - the so-called Ernst equation |22j, [23 
(ft/) Af = (V/) 2 , (1) 

2 The construction of solutions to the Ernst equation by means of Racklurid transformations 
belongs to the p owe rful analytic methods developed by several authors 110 -120 . For a detailed 



introduction see [El 



2 



dp 2 p dp dC, 2 ' \dp' d( 

where the Ernst potential / is given by 

f = e 2U +ib with = a iP , b p = a^. (2) 

The remaining function k can be calculated from the Ernst potential / by a line 
integration: 

^ = 2U P U C + \e~ w b. p b^. 
p I 

To obtain the boundary conditions for differentially rotating disks of dust, one has 
to consider the field equations for an energy-momentum-tensor 

T lk = eu l u k = cr p (p)e u ~ k S(C)u t u k , 

where e and a p stand for the energy-density and the invariant (proper) surface 
mass-density, respectively, S is the usual Dirac delta-distribution, and u 1 denotes 
the four-velocity of the dust material^]. 

Integration of the corresponding field equations from the lower to the upper side of 
the disk (with coordinate radius po) yields the conditions (see (24), pp. 81-83) 

2na p = e u - k (U, c + ^Q) (3) 

e 4U Q 2 +Q(e 4U ) x + (b iP ) 2 = Q (4) 

for C = 0+, < p < po and 

Q = -pe- 4U [b p b x + (e 2U ), p (e 2U ) x }. (5) 

Note that boundary condition ([|) for the Ernst potential / does not involve the 
surface mass-density a p . This condition comes from the nature of the material the 
disk is made of. Therefore, equation ([|) will be referred to as the dust- condition. 
Instead of prescribing the proper surface mass-density a p [which leads to the bound- 
ary condition (|^)] one can alternatively assume a given angular velocity f2 = fl(p) = 
u v ju l of the disk which results in the boundary condition (£ = + , < p < po): 

n = —®—. (6) 

a-X ~ a y 

The following requirements due to symmetry conditions and asymptotical flatness 
complete the set of boundary conditions: 

• Regularity at the rotation axis is guaranteed by 

|(o,C) = o. 

• At infinity asymptotical flatness is realized by U — > and a — > 0. For the po- 
tential b this has the consequence b — > b^ = const. Without loss of generality 
this constant can be set to 0, i.e. / — ► 1 at infinity. 

• Finally, reflectional symmetry with respect to the plane £ = is assumed, i.e. 
f(p, — C) = f{p, C) (with a bar denoting complex conjugation). 



3 u 8 has only ip- and t- components. 
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1.2 Solutions of the Backhand type 



For a given integer q > 1, a set {Yi, . . . , Y q } — {Y,},j[] of complex parameters, and 
a real analytic function g defined on the interval [0, 1], the following expression 



f = fo 



with (a bar denotes complex conjugation) 
l 



■ fa = exp 



(-l) q g{x 2 )dx 
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Zd = y/(ix - C/P0) 2 + (P/Po) 2 , 3?(Z D ) < 



X v = 



Y v -m 



z = -( P +ic), a: a, = 1 



7^ 



7„ = exp A !/ (F„ + iz) 



(-\y g (x 2 )dx 
(ix -Y V )Z D 



a* a„ = 1 



satisfies the Ernst equation. With the additional requirement that for each pa- 
rameter Y v there is also a parameter Y M with Y, = — Y M , reflectional symmetry, 
— = /(PjOj i s ensured^. Moreover, the parameters Y„ are assumed to lie 
outside the imaginary interval [— i, i]. 

The above Ernst potential / = f(p/po, C/Po; {^}<? i <?) is obtained by a Backlund 
transformation applied to the real seed solution / , sec [pH , On the other hand, 
as demonstrated in appendix A, it can be constructed from the hyperelliptic solu- 
tions by a suitable limiting process (see also Q). The particular ansatz chosen for 
the seed solution / guarantees a resulting Ernst potential which corresponds to a 
disk-like source of the gravitational field (see also section 1.2 of |Q). 
Furthermore, / does not possess singularities at (p, C) = Po(|3[Yy]|, — 5J[Yy]). This 
is due to the fact that a v X v is a function of A 2 ,, and this means that / does not be- 
have like a square root function near the critical points (p, £) = pod^Y,,] |, — 5ft[Y„]), 
but rather like a rational function. Now, in the whole area of physically interesting 

4 In the following, the notation {Yi, . . . , Y q } will be abbreviated by {Y^},j. 

5 Hence, the set {iY^},j consists of real parameters and/or pairs of complex conjugate parame- 



ters. 
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solutions that will be treated in the subsequent sections, each zero of the denom- 
inator is cancelled by a corresponding zero of the numerator in (^) such that the 
resulting gravitational field is regular outside the disk. 



The real function g that enters the Ernst potential is assumed to be analytic on 
[0, 1] in order to guarantee an analytic behaviour of the angular velocity Q for all 
p G [0,po]- Moreover, the additional requirement 

.9(1) = 

leads to a surface mass density o~ p of the form 



ap(p) = a 4>p[(p/ 'p ) 2 }\/l - {p/po) 2 [with tp p analytic in [0, 1], ip p (0) = 1] (8) 



and therefore ensures that o~ v vanishes at the rim of the disk. 



In this article the question as to whether the above expression for the Ernst potential 
is sufficiently general to approximate arbitrary differentially rotating disks of dust is 
investigated. Of particular interest is a rapidly converging method to perform this 
approximation. To this end, the set {^jg of complex parameters will be translated 
into an analytic function 

i : [0, 1] -» M. 

Thus the Ernst potential will depend on two real analytic functions defined on [0, 1]: 

/ = f(p/po,t/po;€;g), 

which eventually proves to be sufficient to satisfy both the dust condition (^) and the 
boundary condition (|^) [or alternatively @] . The rapid and accurate approximation 
can be realized since both g and £ are analytic on [0, 1] and thus permit elegant 
expansions in terms of Chebyshev polynomials. 



2 Generalization of the Backlund type solutions 
by a limiting process 

As demonstrated in JlJ for the Backlund type solutions with q = 1, the dust condi- 
tion (Q) can be satisfied by an appropriate choice of the function g if the complex 
parameters Y v are prescribed. To fulfil a second boundary condition, (0) or (^), the 
set {y } q of these parameters has to be translated into a real analytic function £. 
To this end, consider the following equalities for the above solutions / = f{{Y u } q ; g)f\ 
which are proved in appendix B: 

/[{y 1) ...,y,_ a ,y,_i s yj ;fl ] = f[{Y u ...,Y q _ 2 } ]9 ] (9) 

if y 9 _! = —Y 9 e k 

f[{Y u ...,Y q ^,Y q ^,Y q };g} = f[{Y u . . . , Y^}; g] (10) 

if y,_! = Y q 

lim f[{Y 1 ,...,Y q _ 1 ,it};g] = f[{Y u ...,Y q - 1 };g] (11) 

if t G R 

lim f[{Yi,...,Y q - 2 ,Y q -i,Y q };g] = f[{Y lt . . . , Y q ^}- g] (12) 

if y 9 -i = -Y q . 

6 In the following the Ernst potentials / given by (Q) are considered as complex functions 
depending on the set of complex parameters and on g. 
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In order to find an approximation scheme, the desired function £ = (,{{Y u } q ) is 
supposed to be invariant under the modifications (pfjl^) of the set {y„} q that do 
not effect the Ernst potential. This property will be necessary to solve the boundary 
conditions uniquely. 

It is realized by the real analytic function 



ax 2 ;{Y u } q ) = -ln 



1 :y 



x e [-f,f] 



(13) 



which can be proved by considering that for each parameter Y v there is also a 
parameter Y^ with Y v = —Y^ , and that, moreover, the parameters Y v do not lie 
on the imaginary interval [— i, i]. 



The set X of all functions £ = t;(x 2 ; {Y u } q ), q € N, which are defined by ( |l3|) forms 
a dense subset of the set A of all real analytic functions on [0, 1]. Now, for a given 
function g, each £ G X is mapped by (0) onto a uniquely defined Ernst potential 



%:X^£, %(0 = f({Y„} q ;g), 



(14) 



where the set {Y v } q results from £ by (|l3|). 

In the following, it is assumed that this mapping <I> ff can be extended to form a 
continuous function defined on „4.^| Then, given the two real functions g and £, 
defined and analytic on the interval [0, 1], the Ernst potential 

/(£<?)= lim f({Y^} q ;g) 

q — >oo 

exists and is independent of the particular choice of the sequence {{fJ 9 ^ } q}%L qo 
which serves to represent £ by 



— lim In 

X q^oo 



n 



iY,} 



(?) 



= 1 1 1 V 



(?) 



for x £ [-1, 1]. 



This provides the groundwork for the approximation scheme that will be devel- 
oped in the next section. The treatment additionally assumes that the boundary 
conditions @ and (f|) [or (Q) and (||)] interpreted as functions of g and £ are invert- 
ible. The accurate and rapid convergence of the numerical methods justifies this 
assumption although a rigorous proof cannot be given. 



3 An approximation scheme for arbitrary differen- 
tially rotating disks of dust 

It is now possible to attack general boundary value problems for differentially ro- 
tating disks of dust. With the above generalized solutions / = /(£; g) the boundary 
conditions [see formulas (|3||^, ^)] become a problem of inversion to determine g and 
£ from a p or Q: 

(A) S(g;0 = {e u - k [U x + ±Q]/[<Wl - (p/po) 2 ]}(£; g) = 2tt0 p or 

(A>) O( 5 ;f)={g/[O(0)(a, c -aQ)]}fc S )=fl/n(0) = n* (15) 

(B) D(g;Q = {p 2 [Q 2 e w + Q (e* u ) x + (b p ) 2 ] } (£; g) = 0, g(l) = 

7 Here, £ denotes the set of all Ernst potentials corresponding to disk-like sources. 
8 The mathematical aspects of this assumption will be discussed in section 5. 
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This inversion problem is tackled in the following manner: 

1. The only way to treat the complicated system ( |l5| ) numerically seems to be 
by restricting it to a finite, discretized version and solving this by means of a 
Newton-Raphson method. 

2. For this method, a good initial guess for the solution is needed. As shown in 
appendix C.l, there exists a representation of the functions g and £ in terms of 
dp or in the Newtonian regime e <C 1 where e = M 2 /J and the gravitational 
mass M and the total angular momentum J are given by 

M = 2 j (T ab - ^Tg ab )n a fdV , J = - / T ab n a V b dV , T ab = g ab T ab . (16) 
s s 

(S is the spacelike hypersurface t = constant with the unit future-pointing 
normal vector n a ; the Killingvectors £° and rj a correspond to stationarity and 
axisymmetry, respectively.) 

3. This motivates the following finite version which results from expansions of 
dill ) in terms of Chebyshev-polynomials Tj(r) = cos[j arccos(T)]: 

Fj{v k ) = (l<j,k<Ni + N 3 -l): 

■ Fj = Dj (l<j<Nt- 1), F Nl = e(g m ; £») - e, 

i^+j-i = - 2tt^ or = Oj - n; (2 < j < N 2 ), 

v k =g k+1 (1 < Jfe < JVi - 1), ujv 1+ fc_i=& (l<fc<iV 2 ) 

ATi iVi 

• ff(^ 2 ) « E ft^-r(2x 2 - 1) - i 9l , ,g(l) = => 9l = -2 £ # 

3=1 " 3=2 

•£(* 2 )~ E^_i(2ar 2 -l)-|£ 1 
3=1 
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3=1 

V> P (0) = 1 => ^ = 2 § (-l)tyi + 2 
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•^(a: 2 =p 2 /p 2 ;<?;£)« E ^( 5m ;^)T :? - 1 (2x 2 - 1) - i£>i( ffm ;^) 

3 = 1 
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[The function s(g m ]£,n) = M /J is determined using (|16|) for the above func- 
tions g and £.] 

4. For the above system, the boundary values are assumed to be given in the 
form of the ipk's or f2j£'s (k — 2, ...,A^)- Moreover, some e <C 1 has to 
be prescribed. Then, good initial v^s come from the Newtonian expansion. 
The Newton- Raphson method improves the v^s and yields a very accurate 
solution to ( |l5| ) for the chosen small e. Now, this solution serves as the initial 
estimate for the v^s belonging to a marginally increased value for e. Again, 
the Newton-Raphson method improves the solution, and one continues in this 
manner until this procedure ceases to converge. This occurs for some finite 
value £o: at the latest for e = 1. A further discussion of this limit is given 
below. 

5. A rather technical detail is the retranslation of the £j into a set {3^} 9 which 



then gives a satisfactory approximation of £ in terms of (13). There are many 
ways to do this. Here, the following one has been chosen. 
One rewrites equation ( |l3| ) in the equivalent form 

exp[^(, 2 ;{y;} g )]=ni|i T | = ^ with P q ( X )=±b vX ". 



The coefficients b v of the polynomial P q can be determined by evaluating the 
left hand side at q arbitrary different points x p S [0, 1] [] and solving the 
following linear system: 



exp [ Xll £(a£; {Y v } q )\ V b u x^ = V b„(-x^y 



The zeros of P q determine the Y v . 

The above scheme has been performed for many different prescribed surface mass 
densities and angular velocities. This provides strong evidence for the conjecture 
that, in this manner, all Newtonian disks can be extended into the relativistic 
regime. It has been found that the value for £rj, the limiting parameter for the 
convergence of this scheme, depends on the chosen profile for ip p (or equivalently 
for fi*). It is illustrated in appendix C.2, how the Ernst potential always tends to 
the extreme Kerr solution |25| as e — ► 1. This supports a conjecture by Bardeen 
and Wagoner ^6|. But eq = 1 does not hold for all given surface mass densities. 
Even in the Newtonian regime there are surface mass densities for which a realistic 
physical disk cannot be found since the corresponding angular velocity would be- 
come imaginary. If one chooses a profile for a p not very different from these, then 
the Newtonian limit still might exist, but some Eq < 1 turns up, beyond which the 
method does not converge. In the case of prescribed angular velocity, the situation 
is similar. Here, for any sequence / = f(g £ ; £ e ) the angular velocity O* tends for all 
x 1 € [0, 1] to 1 as e — * 1. So, each nonuniform rotation law will lead to some £o < 1 
(see section 4 for examples). 



The above expansions in terms of Chebyshev-polynomials allow a very accurate 
representation with only a small number of coefficients. However, the retranslation 
of £ (see the above point 5) leads to functions that are not especially well suited for 
an approximation. In particular, if the boundary condition ip p is chosen to be close 
to those for which there is no Newtonian disk, then the accuracy cannot be driven 

9 Here, zeros of Chebyshev-polynomials have been used. 
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particularly high by the computer program used, although the method in priniciple 
allows arbitrary approximation (see section 4.2). 

For ^ p 's sufficiently far away from those critical ones, the accuracy obtained was 
very high. By choosing appropriate values for Ni and N2 one can always achieve 
extremly good agreement with the dust condition (Q) (12 digits and beyond) which 
ensures a realistic physical interpretation of the solution. The accuracy to which 
the second boundary condition, (^) or ([j|), can be satisfied, depends on the pa- 
rameter e. It is usually around 8 digits in the weak relativistic regime, and falls 
as e increases, but is still around 4 digits as e tends to eo- These values arose for 
Ni = 30, N2 = 12, and typical ipp's (like ip p 's depending linearly on a; 2 ) and 0*'s 
(e. g. the realistic one considered in section 4.1). The number q of the parameters 
Y v by which £ is represented, was chosen to be between 20 and 30 (independently 
of N 2 )- 



What remains to be discussed is the regularity of the Ernst potentials that were 
obtained. For a few of the solutions, the functions e 2U and b were plotted over the 
coordinates p and £. Moreover, the agreement of the alternative representations of 
M and </, as given by the behaviour of the Ernst potential at infinity 

U = -— +0(r- 2 ), b= -2J^ + 0(r- 3 ), (r = , C = r cos 9) 



with the results from formulas (16) yields good confirmation of the regularity. This 



agreement was checked for all solutions that were calculated. 



4 Representative examples 

From the numerous solutions obtained, three particular sets of differentially rotat- 
ing disks are discussed in more detail. The first one is an example of disks revolving 
with a realistic rotation law. The second set illustrates the break down of the nu- 
merical method for a specially prescribed surface mass density a p at some £0 < 1. 
On the other hand it is demonstrated that, for the same o~ p , regular solutions can 
be found in the highly relativistic regime. Finally, the third example concerns the 
occurence of a second ergoregion for a particular series of disks and, moreover, the 
gradual merging of the two spatially separated ergoregions as e increases. 
The deviations between the boundary values obtained for particular numerical solu- 
tions and the given boundary conditions are listed in tables. The quantities A^j, Aq, 
and A (j therein are defined by 

A D = max \D oht (x 2 ;g;£)\ 
z 2 e[o,i] 



max 

a; 2 S[0,l] 



^obt(* 2 )- fi giv(* 2 ) 



max 

z 2 e[o,i] 



where the indices 'obt' and 'giv' refer to obtained and given quantities, respectively. 
Moreover, by letters (a), . . . ,(e), special examples are marked, for which illustrative 
graphs have been made. Here, curves drawn in the same line style belong to the 
same solution. The graphs show the dimensionless quantities poo~ p and poQ as well 
as g and £ plotted against the normalized radial coordinate p/po and x, respectively. 
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4.1 Disks possessing a realistic rotation law 



As motivated by observations in astrophysics the rotation law of a galaxy is often 
modelled by an equation of the form (see [^7| ) 



fi(p) 



n(o) 



(17) 



Here, the parameter p\ varies for different galaxies. In the following series of solu- 
tions illustrated in figure 1, p\ = 0.7/?o has been chosen. As described in section 3, 
there is a limiting parameter Sq « 0.935, for which the numerical method ceases to 
converge. 
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Figure 1: Disks possessing the rotation law p7| ) with p\ = 0.7po 
(JVi = 30, N 2 = 12). 



4.2 Disks with a critical surface mass density 

For the following sequence of solutions, a surface mass density of the form 

o2 



<r P {p) = o-o(l-3^2+/3^j 
Po Po 



1- 



(18) 
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has been assumed. 
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Figure 



2: Disks possessing the surface mass density jl"8| ) with (3 = 6 
(Ni = 30, N 2 = 12). 



It turns out that for (3 > /3/v « 7 no Newtonian disks with a real angular velocity 
can be found. On the other hand, for (3 = 5.5, all relativistic solutions for < e < 1 
exist. The table and graphs of figure 2 refer to the case (3 = 6. Starting here from 
the Newtonian solution, one soon recognizes a first limiting parameter e$ w 0.60 for 
which the method breaks down. However, by coming from solutions with (3 = 5.5 
and e close to 1, it is possible to create highly relativistic solutions with (3 = 6. In 
fact, there is another limiting parameter, e\ w 0.97, above which the solutions with 
(3 = 6 exist once again. Due to the nearness to the critical surface mass density (for 
(3 = Pn), the accuracy obtained for the boundary condition (H) is not very high. 



4.3 Disks possessing spatially separated ergoregions 

The particular set of disks depicted in figure 3 demonstrates the occurence of a 
second ergoregionQ These solutions do not satisfy a specially prescribed bound- 
ary condition (|^) or (^), but have been constructed in the following manner as 
intermediate solutions. 

10 An ergoregion is a portion of the (p, £)-space within which the function e 2u is negative. 
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Figure 3: Example for a series of disks possessing spatially separated ergoregions. 
In the uppermost picture, the rims of the ergoregions in the (p/po, C/po)-space are 
to be seen (Nx = 40, N 2 = 9). 
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If one investigates solutions with surface mass densities similar to those of (jig), 
one recognizes two minima for e 2U (taken as a function of p, < p < po? C = 0), 
say at p a and pb > p a - Now, for a particular choice of a p it is possible to get 
e 2U (p a ) > and e 2U (pb) < 0, whilst by another choice one can achieve e 2U (p a ) < 
and e 2U (pb) > 0. This makes clear, that disks with spatially separated ergoregions 
can be constructed by interpolating between these solutions. 

For the chosen example, there is only a narrow interval (£ a ,e&) for which the two 
separated ergoregions occur. As can be seen from figure 3, after creation of the 
second ergoregion at e a ~ 0.8403, both ergoregions grow as e increases. Eventually, 
at Eb ~ 0.8415, the ergoregions merge into one ergoregion. 



5 Discussion of mathematical aspects 



As already mentioned in section 2, the assumption that the function $ s introduced 
in ( |l4| ) can be extended to form a continuous mapping defined on A, lies at the heart 
of the above numerical methods. Although this assumption seems to be intuitive, 
it is not trivial. Consider the following example: 
For any analytic function ifi : [0, 1] — > R one finds the equality^} 



lim \ In 

a — >oo < J 



1 



1 



From this it follows that 

,i 1 q 

2 / <f>(t) dt = - lim > In 

Jo x rri 



ip(t) dt. 



q + x(f>{v/q) 



with ip(t) = ±xcj)(t). 



Hence, the function £,(x 2 ) 
y(«) = q 



4>{v/q) 



with 



q-x<j>{v/q) 

2 can be represented by any sequence of the form 
4>{t)dt = 1. 



Since these sequences might be quite different from each other, it is rather surprising 
that all of them approximate the same Ernst potential given by (0) . But this follows 
from the above assumption. 

This already indicates the difficulties which are connected with a rigorous proof of 
this assumption because the Ernst potential is only given in terms of the set {^jq 
and not directly in terms of £. 



A further conjecture is strongly confirmed by extensive numerical investigations: 



For the hyperelliptic class of solutions represented by in appendix 
A, the functions £ and g are given by 



5(z 2 ) 



i 

2x 



In 



n 



sign Y[X V \ A g {x 2 )h{x 2 ) , 



A g (x 2 ) 



Y[(ix-X v )(ix-X u ), A g (x 2 )>0. 



llr To verify this formula one simply expands the logarithms in the form ln(l + e) = e + 0(e 2 ) 
and notes that the resulting sum tends to the Riemann integral of the right hand side. 
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In particular, in this formulation, the solution for the Neugebauer-Meinel- 
disk |J |^/ assumes the form f = /(£;<?) where 

1 x 2 - dMx + C 2 (n) 
2x 11 x 2 + d(ji)x + C 2 (m) ' 

d(//) = V2[1 + C 2 (aO], C 2 (/i) = - v /TT7^, 

A 4 

arsinh[/i(l — x 2 )] , 

and i/ie parameter ji, < < fio = 4.62966184.., is related to the 
angular velocity by 

fi = 2n 2 p 2 e- 2Va , V a = U(p = 0,( = 0). 

As already mentioned, a direct proof of the above assumptions promises to be 
very complicated. But there might be an alternative proof which relies on relating 
a general solution of the Ernst equation to the solution of a so-called Riemann- 
Hilbert problem, see (ll| |l], In this treatment, an appropriately introduced 

matrix function, from which the Ernst potential can be extracted, is supposed to 
be regular on a two-sheeted Riemann surface of genus zero except for some given 
curve, where it possesses a well-defined jump behaviour. The freedom of two jump 
functions defined on this curve corresponds to the freedom to choose £ and g. Now, 
if one succeeds in finding a particular formulation of a Ricmann-Hilbert problem 
in which £ and g are involved, then the final solution for / proves to depend only 
on £ (and g) and not on a particular global representation in terms of {Y^} g . This 
deserves further investigation. 



a- 2 ) = 

g(x 2 ) = 



There is very strong numerical evidence for the validity of both assumptions. For 
various functions £ (and functions g), different representations have been seen 

to approximate the same Ernst potential. In particular, the approximation of the 
Neugebauer-Meinel-solution in terms of Backlund solutions was carried out to give 
an agreement up to the 12th digit with the hyperelliptic solution, which confirms 
both assumptions. 



A The transition from the hyperelliptic solutions 
to the Backlund type solutions 

In this section the Backlund type solutions are derived from the hyperelliptic class. 
The latter is assumed to be given in the form represented in for an even integer 
P > 2: 



V =1 x v 



XPdX 
V(Xj 



(19) 



■v(x) = 


p 

(X + iz)(X - iz) J[{X X V ){X - X v ) 







po 



12 The parameters K v , the upper integration limits IfM , and the integration variable K have 
to be replaced by their 'normalized' values X u = K u /po, X^"> = K^/po, and X = K/po, 
respectively. 
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1 



1 

f (ixVh{x 2 ) dx m 
/ ^ — ) — '- — , < j < p, h : [0, 1] -> K, analytic, 
7 AD 



as defined in (Q) 

The set {L¥„} p consists of arbitrary real parameters and/or pairs of complex conju- 
gate parameters (in order to guarantee reflectional symmetry). The (^-dependent) 
values for the X^ v ' as well as the integration paths on a two-sheeted Riemann sur- 
face result from the Jacobian inversion problem (|2(]) . 

The transition to the Backlund type solutions ([?]) can be obtained in the limit e — > 
by the following assumptions: 

• p = 2q 

• X 2v -x — Y v + s(3 v , X 2v = Y V (l<v< q), {(i v } q arbitrary 
. g(x 2 ) = (-iyh{x 2 )A(ix), A(X) = f[ (X - Y V ){X - Y v ) . 

v=l 

To this end, the above expression for / is rewritten in the equivalent form: 

'° f(2 ^ 1) A{X)dX f x(2 " ] A(X)dx\ f {-iy g {x 2 )dx 



f = cxp 



^ f A(X)dX f A{X)dX \ r 

h\J**>-i ~nxV + J X2 „ v(x) ) J 



Z D 



The Jacobian inversion problem ( |20| ) reads as follows in a similarly rewritten form 
(1<M<<7): 



E 



x(2v 15 A(x)dx r xi2v) A(x)dx \ _ r {-\y g {x 2 )dx 

x ^_ x V(X)(X-YJ + J Xav V(X)(X-Y^ j - J (ix-Y^Zj, 



1 



^ / r x(2v 1] A{x)dx t t x(2v) A{x)dx \ _ r {-iy g { x 2 )dx 

'x*-! V(X){X-Y^) + Jx 2v V(X){X-Yp) j ~J (ix-Y^Zo 



u=l 



-1 



Furthermore 

A(X)dX ; f x '"' A(X)dX 



X-2 



. t V{X)(X-Y) J X2u V(X)(X-Y) 
x ^~ l A(X)dX rx^-v A(X)dX 



X2v V(X)(X-Y) J x w V(X)(X-Y) 

with X( 2l/ ) now lying in the other sheet of the Riemann surface. 
In the limit e — > 0, one obtains 



lim 



'X'"- 1 A(X)dX _ f ±inS llv /[X ll (Y ll + iz)\ for Y = Y /i 
X2v V(X)(X-Y) " o for Y = Y U 
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with <5^„ being the usual Kronecker symbol and A M as defined in (Q). 
The second term amounts to 



lim 



* t2 - I) A(X)dX f^- 1 ' dX 



oJxW V[X){X-Y) J x w (X-Y)y/(X + iz)(X-iz) 



In 



nXjX^-V) - X(Y)] [X(X^) + X(Y)] \ 
X{Y)(Y + iz) ^ \ [A(X(^-i)) + X(Y)} [A(X( 2 -)) - A(F)] J 

where for evaluation of the second integral the substitution 



has been used. 

Hence, the Jacobian inversion problem reads as follows in the limit e — > 0: 

* [X(X^-V)-X,} [X(X^) + X,] 
y i [A(X(^-i)) + A M ][A(X(2-))-A M ] 7 " 1 ; 

« [a(x( 2 - 1 ))-a;] [a(x( 2 -)) + a;] = _ 

• 11 [A(X(^-D)+A*][A(X(^))-A*] 7 " 1 j 

and in an analogous manner 

. _ . A [A(x( 2 - i )) + i] [x(x^y)-i] 

M [A^t 2 "- 1 )) - 1] [A(X( 2 -)) + 1] ( ' 

[with 7^ , A* and f as defined in (Q)] . 

Instead of evaluating the quantities A(X^- 1 ), (1 < v < 2q), the coefficients b v and 
Cjy (1 < v < q) of the polynomial 

i 

P(X) = l[[X- A^ 2 "- 1 ))] [X + X(X^)} 

(24) 

= A^ + A^b^A^+^^A 2 '^ 2 
are determined. Since 

p{K) „ p^D p(-i) r . 



P(-A M ) ,M ' P(-AJ) P(l 

the following system of linear equations for the quantities b u , c v , P(l), and P(— 1) 
emerges: 
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E [ftva'fA-J^-^cCA*) 2 "- 2 ] - -(A;) 2 ' 
v=\ 

£>,-c) + P(-l) = l 



(26) 



5>, + c,)-P(l) = -l, 



with a M and a* as defined in (Q). 

Finally, if the solution of this linear system for P(±l) is expresseed by means of 
Cramer's rule, the desired form (0) of the Backlund type is obtained. 

B Invariance properties of the Ernst potential 

For the proof of the properties the Ernst potential (0) is reformulated by 

D(-l;{Y„} q ;g) 



f({Yv} g ;g) = fo- 



D{l;{Y v } q ;g) 



(27) 



with 



■D(\-{Y v } q ;g) 



a 2 



(aixi) 
{0,2X2) 



fl 2g+l {o-2q+lX2q+\j 



{a\x\ ) 1 X\ 
(a 2 xl~ 1 ) 1 x 2 



(a 2q +ix q 2 +1 ) 1 x 2q +i 



"2q+l 



■0,1 = A , a,2v = ol v \ v , a 2lJ+ i = a* A* , 

■Xl=X 2 , X 2l , = A 2 , X2u+1 = {K) 2 - 

The above expression for D(X; {Y u } q ;g) is a Vandermonde-like determinant. These 
determinants have been studied in detail by Steudel, Meinel and Neugebauer pOfl . 
By their reduction formula [(8) of pQ] ], D assumes the form: 

D(\;{Y„} q ;g) 

= V qtq +i(a r ; b r \ x r ) [with b r = 1 for r = 1 . . . (2q + 1)] 



— ( Il ar 0') ) ^[^(l), ■ ■ ■ ) x r(q)] [x r (g+l) , ■ ■ ■ , 2V(2g+l)] 



where 



• the sum runs over all permutations P — [r(l), . . . , r(2g+l)] of (1,2,..., 2g+l) 
with r{k) < r(j) for k < j < q as well as for q < k < j 
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£p = 



-1 for P even 
-1 for P odd 



• the Vandermonde determinants are given by 

Vn[xi,...,xn]= Yl(xk-Xj). 

k>j 

In this formulation the following properties can be proved: 
(A) If x 2q +i = x 2q then 

D(X;{Y v } g ;g) 



= {-l) q {a 2q - a 2q +i) 



2q-l 
3=1 



{Y v } q -i;—g) 



(B) If x 2q = 1 + ne + 0(e 2 ), x 2q+1 = 1 - Ke + 0(e 2 ), and (a 2q a 2q+1 ) = 1 + 0(e), 
then 



D(=F1;{F*} 9 ;0) 

2?-l 



i=2 



(a 2g + a 2q+1 ± 2)£>(Tl; {^}<z-i; <?) + 0(e 2 ). 



With (A) the equalities @ and ( |To| ) can be derived whilst (B) serves to confirm 
(0) and @. In order to prove (A) consider the following groups of permutations 
separately: 

Pi : r(q - 1) = 2g, r(g) = 2q + 1 

P 2 : r(2q) = 2q, r{2q + 1) = 2q + 1 

P 3 : r(gr) = 2g, r(2q + 1) = 2g + 1 

P 4 : r(?) = 2g + 1, r(2 ff + 1) = 2g 

For X2 9 +i = x 2q , all terms belonging to Pi and P 2 vanish while all terms belonging to 
P3 and P4 possess a common factor, [a 2q Ti^'Li 1 ( x 2q — and [a 2(J+ i n^=i 1 (^z? — 
Xj)], respectively. After reordering (from which the factor (— l) q results), (A) is 
easily obtained. 

The proof for (B) works similarly. Now, eight groups of permutations have to be 
considered separately: 



Pla 


:r(l) 


= 1, r(? - 1) = 2q, r(q)=2q+l 


Plb 


r(q 4 


- 1) = 1, r(2q) = 2q, r(2q + 1) = 2q + 1 


p2a 


: r(q) 


= 2g, r(g + 1) = 1, r(2q + 1) = 2q + 1 


P 2b 


r(l) 


= 1, r(g) = 2g, r(2g + l) = 2^+1 


Pza 


: r(q) 


= 2g + 1, r(q + 1) = 1, r(2<? + 1) = 2g 


Pzb 


r(l) 


= 1, r(g) = 2g + 1, r(2g + 1) = 2q 


P4a 


: r(<? - 


-l) = 2q, r(q) = 2q + l, r(q + 1) = 1 


Plb 


r(l) 


= 1, r(2g) = 2g, r(2q + 1) = 2q + 1 
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All terms of permutations with a coinciding first index can be combined to givep^l : 

{Pia, P2b} => a 2q F + <D{e 2 ) 
{P 3a ,P 3 b}^a 2q+1 F + 0{e 2 ) 
{P 4a ,P4b}=^±2F + 0(e 2 ) 
with 

2q-l 

F = {-\y +1 ne J] (1 - x j )D(±l; {Y v } q _r, -g). 

3=2 

C Newtonian and ultrarelativistic limits 

C.l The Newtonian limit 

In the limit of small functions g and £, i. e. 

g(x 2 ) = e g g (x 2 ) + 0(e 2 g ), Z(x 2 ) = e^o^ 2 ) + 0(e 2 ), 
the Ernst potential / = /(£;g) as introduced in section 2 is given by 



l l 

r f;„\„ (^2\t f„2 



/(£;<?) = 

1 

f g {x 2 )dx f (ix)g (x 2 )Zo(x 2 )dx 2 2 

I- E a J — — ie g e^ J — + 0{e g ) + 0{e g e 6 ). 



(28) 

In this section, the above property will be proved and the functions go and £o will 
be derived as they result from the Newtonian expansion of the boundary conditions. 

C.l.l The Ernst potential for small functions g and £ 

Due to the assumption that the function <& g introduced in ( |l4|) can be extended to 
form a continuous mapping defined on A (see sections 2 and 5), the representation 
of £ in terms of {Vj/} g can be chosen arbitrarily. Here, the following set is 
used: 



q = 4r 

Y<Lv-3 = Z v (l + £^Z U ), Yav-1 




R(z„) ^ o, z v e 

(v = l...r) 



Then, it follows from (|13|) that £(x 2 ) = e^ (x 2 ) + 0(e|) with 



Co(x 2 ) = -4i 



hi (x 2 + Z 2 )(x 2 +X) 



To evaluate the Ernst potential in this limit, the formulation (|21||26|) in appendix 
A is used and the following steps are performed: 

13 Here the requirements a\ = =pl, x% = 1 are necessary. Additionally, for P^ a and P^, the 
constraint a2 q a2q+i = 1 + 0(e) is needed. 
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1. At first, it turns out that in the limit — > the coefficients b v of the polono- 
mial vanish. This can be seen by considering the solution to linear system 



A, 

D 



D 



a 2 



d2q+l 



(a 2 x q 2 x ) 



( a 2 9 +1^2g+l) 1 X 2q+1 



.X 2 



9-1 
C 2q+l 



a 2r] — a v X v , a 2n+ i — a* A* , x 2r] — A 2 , x 2 r l +i — (A*) 2 

is derived from D by replacing the z^-th column by the vector 

{~ x 2 i ■ • ■ 1 ~ X 2q+l\- 

For 1 < v < q , D„ can be expanded in terms of Vandermonde determinants 

V q+ i(x r (i),..., 2V( g+ i)), r(?7) G {2, ...,2?+ 1}, r{r)) < r{y) iox 7) < n . 

In the limit — > 0, any set {av^jjq+i contains at most g different values, 
and therefore all D u vanish. On the other hand, D remains finite (here only 
Vandermonde determinants V q are involved), and hence all &„ tend to zero. 

2. Thus, with any zero A„ of the Polynomial J24]), (— X v ) also becomes a zero as 

— > 0. This set of zeros is ordered in the following way: 

{A(I (1) ), -X(X^), \(Xto-V), -X(X^)} = {X 1 ,-X 1 ..., X g , -A,}, 
Suppose there is a A p different from all zeros : 
A p ^ XiX^-V) = X{X^) and A M ^ -A^ 2 "- 1 )) for all v = 1 . . . q. 

Then, since 7 M ^ — 1 for small g, ( pl| ) cannot be satisfied. 

3. This gives rise to the following ansatz (v — 1 . . . q): 

by which the system (|2l|/[22|) can easily be solved to get the set {k„}2 9 • 



4. Finally, if g(x ) — e g g (x ) + 0(e g ) is considered, then fl28j ) follows from ( |23| ) 
by inserting the values obtained for {X(X^>)} 2q . 

C.1.2 The functions go an d £o as resulting from the boundary conditions 

For any family of Ernst potentials / = /(ff e ;£ e ) describing a sequence of differen- 
tially rotating disks of dust with the parameter e = M 2 /J [M and J as defined in 
@], the following expansion is valid (see pp. 83-89): 

/ = 1 + e 2 ( P , Qe 2 + ib 3 (p, C)e 3 + 0(e 4 ). 
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By comparison with (|2|) one gets 

2 



1 1 

•e2(p,C) = -/ — ^ — , b 3 (p,() = -l ^ ■ 

-l -l 

If the boundary conditions, 

■ a p (p) = a M(p/Po) 2 Wl - (P/Po) 2 e 2 + 0(e 4 ) (with 2 (O) = 1) or 

■ - r> r!i[( /9 / j0 o) 2 ]£ + o(£ 3 ) (with ni(o) = i), 

are given, then it follows from equations (|^-^) that 

■ (ea),c = 4^(70^2^1 _ (p/po) 2 or ( e 2),p = 20^1 P and 
•(63),/) = 2ptt Qi (e 2 ),c- 

By expressing e 2 and 63 in terms of go and Co in these equations, one gets Abelian 
integral equations for £0 and go. Their solutions read as follows: 

go(x 2 ) = -4er (l-2; 2 ) / (sin 2 </>) ip 2 (cos 2 4> + x 2 sin 2 (j>) d(j> 

Jo 

9o{x 2 )£o(x 2 ) = 8<j n (l-x 2 ) (sin 2 tf>) Cli (cos 2 tf> + x 2 sin 2 </>) dej) 

Jo 

[with Oi(x 2 ) = n^x 2 )^ 2 )}- 

Note that only one of the functions "02 and 17 1 can be prescribed since both represent 
different boundary conditions of the same Newtonian potential e 2 . Likewise, the 
constants o~q and Oq depend on each other. Moreover, these constants in terms of 
ip2 and Sli are prescribed by the equation e = M 2 / J . 



C.2 The ultrarelativistic limit 

It is difficult to relate the functions g and £ of an Ernst potential / = /(<?;£) to 
its physical properties like M and J. Nevertheless, if a sequence /(</ e ;£ e ) can be 
extended to arbitrary values e < 1, then, in the limit e — ► 1, the universal solution 
of an extreme Kerr black hole is reached. It is illustrated how this limit results from 
the form (Q) of the Ernst potential. 



If the limit po — > is considered for finite values of r = \/p 2 + ( 2 , then by using 
the formulation (p7f) one gets (with £ = rcos6>): 



f= M_- (-iy g ( x -)dx + o( P 2 ) 



E x r + Po [E 3 cos6 - (-l) q E 2 ] 
E ir + p Q [E 3 cos8 + {-l)iE 2 ] 



o( P 2 ] 



The Ej do not depend on p and £ but on g and £. In particular: 



•-El 



61 (&i2i) 



&2g (b 2 qZ 2q ) 



(hZr 1 ) 1 ^1 



(fe,^- 1 ) 1 Z 2q 



zf- 1 



71-1 
J 2q 
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1 z 1 



z'l 



■ E-2 — 



?2q 



(b2qZ 2q ) 



(b 2q Z. 



2q I 



1 Z 



2q 



7 q 
^2q 



■b 2v - 



tanh 



1 ? {~l) q g{x 2 )dx 

2 / Sx-Y v 



■ Z^v-x — Y v , Z<i v — Y v . 

Clearly, if E\ ^ then limp _>.o / = 1- The Ernst potential passes to an ultrarela- 
tivistic limit if E\ and po tend simultaneously to zero such that0 

po^o 2E 2 po 
exists. Then one gets 



fir 



2S!rjr + E A cos 9 - 1 
2flur + E 4 cos 9 + 1 



The only Ernst potential of this form which is asymptotically flat and regular for 
r > is the extreme Kerr solution. The constant flu is then real and describes the 
'angular velocity of the horizon'. Moreover, J = 1/(40^) = M 2 , and hence e = 1. 
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